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Abstract 

It was recently proposed that the leading singularities of the S-Matrix of AT = 4 
super Yang-Mills theory arise as the residues of a contour integral over a Grassman- 
nian manifold, with space-time locality encoded through residue theorems generalizing 
Cauchy's theorem to more than one variable. We provide a method to identify the 
residue corresponding to any leading singularity, and we carry this out explicitly for all 
leading singularities at tree level and one-loop. We also give several examples at higher 
loops, including all generic two- loop leading singularities and an interesting four- loop 
object. As an example we consider a 12-pt N 4 MHV leading singularity at two loops 
that has a kinematic structure involving double square roots. Our analysis results in 
a simple picture for how the topological structure of loop graphs is reflected in various 
substructures within the Grassmannian. 
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1 Introduction and Review 



A proposal was recently made that all of the leading singularities of M = 4 super Yang-Mills 
theory in the large N limit arise as the residues of a contour integral over a Grassmannian 
manifold [1]. It has been conjectured that these leading singularities may be sufficient to 
determine the perturbative S-Matrix of the theory [2], [3], and this has been confirmed for all 
one- loop amplitudes [I]- [I] and for a few simple examples at higher loops [8]- [10]. Thus it is 
hoped that this strikingly new portrayal of the S-Matrix may be part of a new description 
of scattering, where the extreme simplicity of the S-Matrix itself takes center stage and 
space-time locality is encoded in a complicated way. 

The Grassmannian contour integral was discovered through investigations [H]-[l6j of 
scattering amplitudes and the BCFW Recursion Relations [Ej-[21] in twistor space [25]-[28]. 
inspired in part by the twistor string [29J, but it remains a mysterious new object without 
any clear first-principled derivation. The case for its validity was based on two sources of 
evidence, in addition to the fact that it possess all of the required symmetries [I], [30], [31], 
including dual conformal invariance [32] -[45]. The first piece of evidence was the explicit 
computation of various residues and their subsequent identification among known leading 
singularities pQ, [31]. The second and perhaps more interesting piece of evidence was based 
on an analysis of the residue theorems that follow from generalizations of Cauchy's theorem 
to more than one variable. It was shown in many examples [1] that these residue theorems are 
directly related to space-time locality, as they enforce the cancellation of unphysical poles in 
and the symmetries of tree amplitudes and the Infrared consistency of one-loop amplitudes. 
Some of these residue theorems imply non-trivial relations that do not follow from the one- 
loop IR equations [46] and that were conjectured to follow instead from IR consistency at 
higher loops. 

The purpose of the present paper is to provide a simple picture for how leading singular- 
ities emerge as the residues of the Grassmannian contour integral, which we will refer to as 
£ n ,fe. Our methods allow us to identify a residue of C n ^ corresponding to any given leading 
singularity. We will carry out this procedure explicitly at tree level and at one-loop, and give 
a few illustrative examples at higher loops. Our analysis will be 'kinematical' as opposed 
to 'dynamical' in a sense that will be made clear below, so we will not actually prove that 
every leading singularity is in fact a residue, but we believe our analysis is nevertheless very 
powerful. We find it especially striking that the topological structure of the loop graph cor- 
responding to a given leading singularity is reflected in the structure of the Grassmannian; 
this can be seen already in figure HJ 

Leading singularities and Grassmannian contour integrals are not widely known, so we 
will briefly review both. The computation of scattering amplitudes in terms of their leading 
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Figure 1: We illustrate the way that one-loop leading singularities correspond to certain 
subspaces of the Grassmannian. The object on the left is a one-loop leading singularity, or in 
other words it is the product of four tree amplitudes evaluated on the kinematics determined 
by the quadruple cut of the loop integral. The rectangles in the picture on the right are the 
non-zero entries of the kxn matrix characterizing the Grassmannian; each rectangular block 
shares one row with the block adjacent to it. It should be noted that only GL(k) invariant 
statements about this matrix are physically meaningful. 




singularities is a descendent of generalized unitarity techniques [17] . As will be familiar from 
Feynman diagram computations, loop amplitudes involve various logarithms, dilogarithms, 
and so on that are themselves functions of the kinematical invariants of the scattering process. 
These functions have branch cuts, and one can compute the discontinuities across these 
cuts. Those discontinuities may themselves have branch cuts, and we can compute these 
discontinuities, and so on, until we are left with some pure rational functions (we get many 
different rational functions depending on which branch cuts we use, and which loop order 
we are at). These rational functions are the leading singularities of a scattering amplitude, 
and it has been conjectured [3] that the leading singularities are sufficient information to 
reconstruct the S-Matrices of Af = 4 SYM and M = 8 Supergravity. 

At this point the leading singularity may seem like a rather technical construction, but 
in fact it is a simple and physical object. The reason is that the branch cut of an integral 
(such as a loop integral) is approached when a parameter in the integrand forces the contour 
of integration to encircle a pole. In a local quantum field theory, poles in the loop integrand 
can only come from propagators, so by isolating the discontinuity across a branch cut we are 
forcing the virtual particles in the loops to go on-shell. Leading singularities arise when all of 
the loop integrations are fixed (or 'cut') by the requirement that various intermediate particles 
are on-shell. Thus leading singularities are simply products of tree-level scattering amplitudes 
evaluated with very special kinematical configurations. If the full S-Matrix is determined by 
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leading singularities, then it is determined by the classical scattering amplitudes of the theory 
in the simplest possible way. 

Now let us describe our Grassmannian contour integral. A Grassmannian manifold 
G(k, n) is the space of k dimensional planes in an n dimensional space. A convenient way 
to parameterize the points of G(k,n) is with a k x n matrix C aa , where a — 1, ...,k and 
a = 1, ...,n; the rows of this matrix span a k plane. Note that different C matrices related 
by a GL(k) transformation C aa — > L£Cp a correspond to the same fc-plane, so GL(k) is a 
"gauge symmetry" of our description of the Grassmannian. 

In what follows the parameter n will always correspond to the number of particles in a 
scattering amplitude or leading singularity, and k will represent the total number of negative 
helicity gluons in an all-gluon amplitude (or more generally the R-charge sector), so MHV 
amplitudes [48] correspond to k — 2. The contour integral we will consider is an integral over 
the C matrices with a very special integrand: 

(12-t)(23-(t + l))"..(nl-(t-l)) n^(^-) « 

The factors in the denominator are the determinants of the sequential k x k minors of C, 
explicitly they are 

( mi ...m k ) = e a '- ak C miai ...C mkCtk (2) 

The other piece of the integrand is a product of k superconformal delta functions, and this 
is where the dependence on the kinematic variables of the external particles enters. We 
represent the kinematics with twistor variables W where 

W = (X,M) (3) 

and p, is the Fourier conjugate to the spinor variable A, with = Xa^X. Note that these 
super twistor variables W are in the fundamental representation of the superconformal group 
PSU(2, 2|4). The anti-commuting fj variable is an on-shell superspace coordinate [49] . The 
use of twistor variables for scattering amplitudes has been extensively and pedagogically 
discussed in [14], and on-shell superspace in [3]; we will not review them further here. 

To begin to better understand & let us count the number of integration variables in 
momentum space. To go to momentum space we just Fourier transform with respect to the 
fi a variables, giving 

^^)-J (12 . . .^:.tr+ !» °r:l^:il - 1) f v-w-iK - o.) (4, 
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where the p a are extra spinor variables to be integrated over. We see that after eliminating 
these extra spinors there are 2n delta functions, but 4 of these encode momentum conserva- 
tion. This means that 2n — 4 of the coordinates in the C aa matrix will be fixed by these delta 
functions. Also, some k 2 of the coordinates can be eliminated by fixing the GL(k) gauge 
redundancy of the Grassmannian. All of the remaining (n — k — 2)(k — 2) coordinates are 
free, so L n ^ should be regarded as a contour integral in this many variables. The choice of 
contour or residue can be viewed as providing equations that fix the integration variables, 
but we can perform the contour integral and solve the delta function constraints in whatever 
order we prefer. In [1] we solved the delta function constraints first, and only then performed 
the contour integration, but we will find the opposite order to be more enlightening in what 
follows!. 

Once the contour integration is performed so that we are left with one particular residue, 
the full Grassmannian will be reduced to some 2n — 4 dimensional algebraic subspace param- 
eterized by a highly constrained C aa matrix. As a very concrete example that we will derive 
below, the matrix 

1 c 27 c 28 \ 

c 42 c 43 1 (5) 
c 63 c 64 c 65 1 
c 85 c 86 c 87 1 / 

corresponds to a one-loop leading singularity with a 4-pt MHV amplitude at each of the 
four corners of the 'box' pictured in figure [TJ This is a rather remarkable result, because it 
means that all leading singularities essentially only depend on kinematic invariants through 
2n — 4 special parameters, whereas we might expect them to depend on the n(n — 1) invari- 
ants (ij) and This is especially surprising when we remember that this is an Af = 4 
supersymmetric result, so it holds for all of the various helicity combinations. 

The methods we will develop in the following sections will allow us to pick out the special 
subspaces within the Grassmannian that give rise to any given leading singularity. We will 
show that there is a very simple way to glue together many smaller copies of L n ^ so that they 
sit as subspaces of a larger Grassmannian, where the smaller copies are to be interpreted as 
tree amplitudes (or general leading singularities) at the vertices of a loop diagram that has 
been 'cut' to make a larger leading singularity. 

Once we have identified an appropriate subspace within the larger Grassmannian, there 
still remains the question of whether this subspace can actually be obtained as a residue. We 
show that this is extremely plausible in the appendix. However, our analysis is 'kinematical' 

2 One might worry that there exist contours of integration that are incompatible with the delta function 
constraints. We will never be led to such 'bad' contours, although they are a reasonable motivation for solving 
the delta function constraints before performing the contour integration. 
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as opposed to 'dynamical' because we are not able to actually compute these residues in 
general. A full proof that all leading singularities are residues of C n ^ would require this 
computation, and this is beyond the scope of the present work. 

In the next section we show how leading singularities can be written in twistor space, 
and in particular how they can be computed by 'gluing' together other leading singularities. 
Then in section three we begin by motivating our analysis, and then we proceed to identify 
all tree and one-loop leading singularities. At the end of section three we give some very 
non-trivial higher loop examples, including all generic two-loop leading singularities and a 
four loop object with an interesting topological structure (as a loop graph). Also, to show the 
power of our method we provide an explicit 12-pt N 4 MHV two-loop example whose kinematic 
structure involves square roots of square roots. With section four we conclude and discuss 
future directions. In an appendix we give some details of the computation of the residues 
themselves, including an argument for the existence of the tree and one-loop residues, and 
we give an explicit solution for the NMHV (k = 3) sector. 

2 Leading Singularities in Twistor Space 

Twistor variables are an elegant representation of massless on-shell states, so phase space 
integrals such as 

J d 4 £5(i 2 )M 1 (£)M 2 (-£) (6) 

can be written very simply in twistor space as 

J D 3 WpM 1 (Wp)M 2 (W P ) (7) 

This is an instance of the well-known Penrose transform (25]. It is essentially guaranteed by 
Lorentz invariance and the kinematics of twistor space - in other words, since twistors fully 
parameterize light-like states, what else could an integral over twistor space be but a dLIPS 
integral - but let us derive the result explicitly. 

We begin by recalling that the momentum vector can be written in spinor language as 
the 2x2 matrix 

w*=(£ t-) (8) 

We will use (2, 2) signature to facilitate calculation, but all of the results we will obtain 
can be analytically continued back to the usual (3, 1) Minkowski signature. Now I 2 is the 
determinant of £ ■ a, so we can re-write the phase space integral as 

J dHZ>{l xi l 22 - ^ 2i )M x {i)M 2 {-i) (9) 



5 



and we can do the integral by, say, integrating over £^ to give 
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If we parameterize the remaining integral with 
oo, then the integral becomes 



A a Apa and allow A„ to run from — oo to 



/ 



d 2 XDX P M 1 (X, A P )M 2 (A, X P ) 



(11) 



where DX P = (A dX) is the projective measure on RP 1 . It is easy to go from this spinorial 
representation of the integral to twistor space. If we Fourier-represent the dependence of M 1 
and M 2 on A, we find 

DX P d 2 Xd 2 ^ l d 2 fi2e i[ ^-^' x] M 1 {X^ 1 )M 2 {\fi 2 ) = J D z W P M l {W P )M 2 {W P ) (12) 

with the projective twistor variable Wp = (A, jl). This is the result we wished to obtain. 

We will now make use of the twistor transform in order to represent leading singularities. 
A one-loop leading singularity 




is given in momentum space by 



r 4 

/ H dH^M^tu -£ 2 , ...)M 2 (£ 2 , -£ 3 , ...)M 3 (£ 3 , -£ 4 , ...)M 4 (h, -h, ...) (13) 
J i=i 

where we are including the momentum conserving delta functions in the tree amplitudes 
M.i. The 16 integration variables in the £{ are completely fixed by momentum conservation, 
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which provides 12 constraints, and the condition that if = 0, which provides 4 constraints. 
The ii may in general become complex, and we define the integral in this case by analytic 
continuation^. 

In maximally supersymmetric theories we must also sum over the helicities of the particles 
running in the loop; this is accomplished by integrating over the on-shell superspace variables 
7] or rj [49] ; for extensive discussions and examples of that formalism see [3] . In twistor space 
the one-loop leading singularity of Af = 4 super Yang-Mills turns into the superconformal 
integral 

r 4 

/ J] D 3 ^w t Mi(m, m, ...)M 2 (m, m, ...)M 3 (m, m, ...)a*4(w 4 , Wi, ...) (14) 
J i=i 

where W = (W, rj). This has a simple diagrammatic representation as 




where we have not explicitly indicated the external states of the amplitudes A4i. 

It may seem that we have not made much progress, since we have merely substited twistor 
space integrals for phase space integrals. However, the twistor space integrands will always 
be delta functions, so performing the twistor space integrals will only involve some simple 
linear algebra and book-keeping, making them vastly simpler than their momentum space 
counterparts. This fact is an enormous advantage, and it will allow us to begin to unravel 
the structure of the Grassmannian contour integral C n ^. 

It is straighforward to write higher-loop leading singularities in the same way - beginning 
with some L loop diagram with 4L propagators, we simply replace each propagator with a 
twistor variable VK, and integrate over it. Those familiar with 'Hodges Diagrams' |llj-[T4"] 
may find the picture above familiar, as it is a sort of generalization of those diagrams. In 
fact, our diagrammatic representation of leading singularities is in some sense a realization 
of Hodges' idea of 'twistor quilts' [12] for loop amplitudes. 

3 We will not delve into this issue in detail because it will not be relevant for our analysis, but a more 
precise definition involves re-interpeting the original loop integral as a contour integral around the four l/£f 
poles 
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3 Unraveling £ n ^ 

In pQ it was conjectured that the residues of the multi-dimensional contour integral 

(12 • • • A;) (23 • • • (k + 1) ) • • • (nl • • • (k - 1) ) ^ 

are in one-to-one correspondence with the leading singularities of the S-Matrix of M = 4 
super Yang-Mills theory. This conjecture was based in part on evidence accumulated by 
explicitly computing residues and then identifying them with known leading singularities. 
A proof of this conjecture would require a specification of the residues of L n ^ along with a 
'dictionary' relating them to the leading singularities. In this section we will show how any 
leading singularity can be identified with a residue of £ n ,fc- 

In order to relate leading singularities to residues, we need a way to label them both. 
A leading singularity can be specified by drawing an L loop diagram with 4L propagators. 
When each of these propagators is cut, we will be left with a product of tree amplitudes 
evaluated with very special kinematics. If these tree amplitudes are MHV or anti-MHV, 
then we have a single term, or a 'primitive' leading singularity. Otherwise, we will have a 
sum of terms, and although one can regard this sum itself as a leading singularity, it is the 
individual terms in the sum that are residues of £ n ,fc- So we should proceed to write each 
tree amplitude as a sum of terms via the BCFW recursion relations; choosing any one term 
from each tree amplitude gives a primitive leading singularity. 

This last step in the definition may seem a bit arbitrary, but fortunately it can be given a 
nice interpretation. As originally shown by Britto, Cachazo, and Feng [T7j, each term in the 
BCFW recursion relations can be interpeted as the quadruple cut of a one-loop box (if the 
tree amplitudes at the corners of the box are all MHV or anti-MHV, this is just a one-loop 
leading singularity). This means that wherever we see a non-MHV tree amplitude, we can 
replace it with a sum over quadruple cuts of one-loop boxes. This process expresses an L loop 
object with 4L cut propagators in terms of an L+ 1 loop object with 4L + 4 cut propagators. 
If we repeat the process until it terminates, we will be left with a unique product of MHV and 
anti-MHV tree amplitudes at L + SL loops evaluated on the kinematics specified by cutting 
the 4L + 45L propagators. Thus each and every term in a leading singularity computed at 
L loops is in fact itself a leading singularity at L + 5L loops. 

We must also label the residues of The denominator of L n) k is a product of n 

determinants, so in simple cases it is sufficient to specify on which of these determinants 
we are evaluating the residue (or in other words, which factors in the denominator vanish). 
However, for even moderately large n and k this is inefficient because the residues are highly 
'composite' [I], meaning that not only the determinant factors in the denominator vanish, but 
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also their derivatives, second derivatives, and so on. Furthermore, the equations that follow 
by requiring that these determinants vanish can have a very large multiplicity of solutions, 
so this method of labeling does not specify a unique residue. 

This line of thought suggests a better way of labeling the residues. A residue is given 
by solving a large system of algebraic equations for coordinates on the Grassmannian, so 
it is natural to label the residue by the solution itself. In particular, since points in the 
Grassmannian can be specified by a k x n matrix C aa modulo a GL{k) gauge redundancy, 
it is natural to label residues by specifying the explicit form of C . Naively this sounds like 
it could be very involved, since one might expect complicated algebraic relationships among 
the Grassmannian coordinates. However, we will see that even in very general cases the C 
matrix takes a form that is both simple and transparently connected to the physics. For 
instance, in the case of tree level and one-loop leading singularities we will see that the 
C matrix can be fully specified by stating which of its entries are zero in a particularly 
convenient GL(/c)-gauge. We will also see that the topology of the loop diagram representing 
the leading singularity is beautifully reflected by its corresponding Grassmannian locus. 

Although we will show how to identify a residue of L n ^ corresponding to any leading 
singularity, our analysis will not result in a complete proof that these leading singularities 
are actually given by the residues in question. The deficit is due to our inability to compute 
general composite residues. This one remaining issue is a precise mathematical problem with 
a known answer, but its solution should be physically interesting, as the computation of 
composite residues contains most of the dynamical information of C n ^- 

3.1 A Simple Tree-Level Illustration 

In [H], [15] it was shown that tree level scattering amplitudes in M = 4 super Yang-Mills 
theory become very simple when transformed to twistor space. These twistor transformed 
amplitudes gave way to new expressions for amplitudes in both twistor space and momentum 
space using the so-called 'link representation'. As an example, the 6-pt NMHV amplitude 
can be expressed as a sum of terms of the form 



in the link representation, where we are ignoring an overall sign factor. For our purposes, the 
only thing to notice about this formula is that C52 is being set to zero by a delta function. 

In [TJ we described the contour integral which we conjectured contains all the leading 
singularities in the M = 4 theory as its residues. We first discovered this formula by trying 
to interpret 6(052) not as a delta function but as a contour integral around the pole I/C52. In 




(16) 
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fact one can write 




dcije 



.icijWi-Zj 



C52C36Ci4(Ci2C54 — C14C52) (C14C36 — C16C34) (c^qC^ ~ C32C56) 



1 



(17) 



and observe that £g,3 reduces to the U above on the residue of the pole C52 = 0. The discovery 
of C n ± was motivated by a desire to understand how the locality of the S-Matrix is encoded 
in efficient, on-shell methods such as the BCFW recursion relations, where locality seems 
to be quite obscure. In fact as shown in [lj locality is encoded via the very many residue 
theorems that relate the various residues of £ ni £@. 

However, now that C n ^ is known, we can reverse the historical logic. We know that the 
U above is a term in a 6-pt NMHV tree amplitude, so we could use its explicit form in the 
link representation to determine which residue of ^6,3 it comes from. In what follows we 
will unravel the embedding of leading singularities among the residues of by identifying 
them with (very general) link-representation formulas. In the following three sections we will 
recursively identify as residues all the one-loop and tree-level leading singularities of N = 4 
super Yang-Mills, and then explain how the method generalizes to arbitrary loop order. In 
the appendix we use our method to give an explicit formula for all NMHV (k = 3) residues. 

3.2 All One-Loop Leading Singularities 

Now we will use what we have learned to identify the residues corresponding to all one-loop 
leading singularities. To do this we need only compare the expression for Cn,k with the 



integral 




i=l 



(18) 



This integral can be visualized as the diagram 




r 



J 



4 For a different and very interesting approach to this question see [50] and also [16] . 



10 



where we are integrating over the Wi with % = 1,2,3,4, which correspond to the on-shell 
intermediate propagators in momentum space, and each D has rii — 2 external particles that 
are not explicitly displayed. We have labeled the W ai with an index a; where i = 1,2,3,4 
denotes the particular £ l n . ki to which it belongs, and we have a total number of particles 
N = rii + n 2 + n 3 + n 4 — 8 and number of negative helicities K = ki + k 2 + k 3 + A; 4 — 4. Since 
the four D depend on the W variables only through S^ 4 (c aa W a ), performing the integrals is 
a matter of book-keeping. 

We will choose to only partially fix the GL(ki) redundancy of the matrices C l a . a . so that 



CI 



1 1 




Q,(n-2)i 


o\ 





Cai,U ■ 


• • C ait (n-2)i 








Coti,li ■ 


■ ■ C ai ^ n -2)i 





v° 


C i+l,U ■ 


■ ■ Cj+l,(n-2)i 





(19) 



or in other words, we have fixed the first and last columns of the matrix, which correspond 
to the Wi variables over which we are going to integrate, but we have not fixed the other 
columns and rows. There is a subtlety when treating the anti-MHV 3-pt amplitude, because 
its C matrix has only a single row and therefore we can only fix it to be C = (1, cn { , C12). 
We will return to treat this special case at the end. 

We will choose to use a delta function from D to perform the integral over W{. Naively 
one would expect to simply solve for the Wi, but the twistor variables are projective, so we 
can only conclude that 



Wi 



a; 



~Ci aiWa; 



(20) 



for some non-zero Tj. This new Tj is an arbitrary parameter, so we can use it to fix one of 
the c variables, so e.g. we could set c ii i i = 1 for each i = 1,2,3,4. This follows because we 
can then absorb Tj everywhere else it appears by re-scaling the other variables. However we 
will ignore the Tj for now in order to avoid breaking any symmetries. Wi appears in both D 
and C 1 ^ 1 ; subsituting it into the latter takes 



S^id^Wa^+Wi) - S^Ci^Wa, 



i,aiW ai ) 



(21) 



and now we are done! The one-loop leading singularity corresponds to Cn,k with C matrix 
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fixed to the form 
/ c i,ii ••• 

Ca,li • • • 

Ca,li • • • 
Ca,li • • • 

... 



c 










V 



c l,(n-2)i 
C Q ,(n-2)i 

C Q ,(n-2)i 
C a ,(n-2)i 



















C 2,l 2 
Ca,l 2 

*-'Q!,l2 
















C2,(n-2) 2 
c a,(n-2) 2 

C Q ,(n-2) 2 


















C 3,l 3 

Cq,13 
Cct, I3 

C»,l3 













C 3,(n-2) 3 
C a ,(n-2) 3 

Ca,(n— 2)3 
Ca,(n— 2)3 






C a,l 4 














C4,l 4 
Ca,l 4 

c a,l 4 



Co,(n-2) 4 \ 












c 4,(n-2) 4 
C Q ,(n-2) 4 

C Q ,(n-2) 4 / 



(22) 



where we have eliminated the minus signs in front of the Ci jQ>j variables by a simple redefinition. 
We did not completely fix the GL(ki) gauge redundancies of the D ni k . in order to avoid 
obscuring the structure of this matrix, but in practical computations one would fix these 
redundancies in some way. Also, although we have written the matrix as almost-block- 
diagonal, the diagonal of the matrix plays no special role - we are free to cyclicly permute 
the columns and rows. We should think of this C matrix 'picture' as a specification of the 
linear dependencies among its various columns. 

Let us count the number of free variables in momentum space to show that the contour of 
integration has been completely specified. After fixing the GL(ki) redundancies and choosing 
a particular residue for the D ni k . we are left with 2rii — 4 variables in each D [lj, which would 
be fixed by delta functions were we to transform back to momentum space. This means that 
there are a total of 2N free variables after the individual D contours have been specified. 
However, we saw above that there are four T{ parameters which appear as a consequence 
of the fact that we have integrated over R 4 instead of RP 3 four times; we can use these to 
eliminate four c variables by setting them to 1. If we take Cn^k to momentum space we find 
2N — 4 delta function constraints, which is exactly equal to the number of free variables. 

Before giving some examples let us return to the case where one of the C\ say C 1 , is an 
anti-MHV 3-pt amplitude. Let us fix its C 'matrix' to be 

C= (l,c lll)Cl2 ) (23) 



so that the amplitude becomes 



^3,1 



dc\ 1 , dc 



lli«<-12 r-414 



Cll 1 Ci2 



5^(W 1 + c lll W ll +c l2 W 2 ) 



(24) 
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In accord with our choices above we will use this delta function to integrate over Wi, giving 

m = ni-cn.m, - c 12 W 2 ) -> c ni W u + W 2 (25) 

with an appropriate choice of the free parameter t\ and re-scaling of cnj. Now we have 
completely eliminated C 1 , its only remnant being the Cn 1 parameter. As before, we will 
solve for W2 using a delta function from £ 2 , so the end result is a C matrix for Cn,k of the 
form 



/ Ci,li c 2j i 2 
c Qj i 2 



C 











V 















C2,(n-2) 2 

C Q ,(n-2) 2 

C a ,(n-2) 2 

Ca,(n-2) 2 c 3,l 3 

c a ,i 3 








C3,(n-2) 3 
C Q ,(n-2) 3 



c a,l4 



















c a,l 3 
C a ,l 3 






C Q ,(n-2) 3 
C a ,(n-2) 3 C 4 ,l 4 

c a ,i 4 







^0,14 



C Q ,(n-2) 4 \ 











C4,(n-2) 4 
C Q ,(n-2) 4 

C a ,(n-2) 4 / 



(26) 



Let us now check these very general results with a few examples. If we want to obtain a 
box coefficient (one-loop leading singularity) in the MHV sector, we must make one pair of 
opposite corners MHV and the other pair anti-MHV 3-pt amplitudes. This gives a C matrix 
structure 



C = 



* 



* * 

* * * 



(27) 



where there are still two r parameters to be specified (in other words, we can rescale the two 
rows independently by an arbitrary factor, setting a c parameter in each equal to 1). 

One might wonder what would have happened if we made the two anti-MHV 3-pt ampli- 
tudes adjacent. Physically, this sort of leading singularity must vanish; our results give 



C = 



* * * 
0* 



(28) 



In this case the sub-determinant (1,2) vanishes. If we interpret this as 1/0 it means that 
our result is not well-defined. If we attempt to view £ n 2 as a contour integral evaluated 
on the residue (1,2), then when we return to momentum space we would find an additional 
constraint on the momenta beyond momentum conservation, or in other words we would find 
that this object vanishes for generic momenta. Thus we see that C n ^ "knows" that this is 
not a viable leading singularity. 
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Finally let us consider a much more non-trivial example. In the case N = 8, K = 4 there 
is a single four mass box which corresponds to rii = 4, ki = 2 for all i, or in other words this 
is a box with a 4-pt MHV amplitude at each corner. Eliminating the four extra variables, 
we obtain a matrix structure 

/ * 10000* * \ 
* * * 1 
0***100 

yoooo***iy 

We immediately see that the determinants (1,1 + 1,1 + 2,1 + 3) vanish for I odd but that 
they are non-vanishing for / even. This was precisely the residue found in p] to correspond 
to this particular leading singularity. 



C 



(29) 



A Worked Example 

In the analysis above we saw how one-loop leading singularities correspond to particular C 
matrix structures, or in other words, to particular subspaces of the Grassmannian. However, 
we did not show how one obtains these C matrices from contour integration, and we did not 
work out the resulting residues. We will go through these procedures in detail for the n — 8, 
k = 4 example, and then we will explain how they generalize. 
We would like to fix the GL(4) redundandancy so that 

/ c 2 i 1 x 3 c 27 c 28 \ 

c 4 i c 42 c 43 1 x 5 

c 63 c 64 c 65 1 x 7 

\ xi c 85 c 86 c 87 1 

However, there is a non-trivial Jacobian that arises when we fix the GL(4) redundancy in 
this way. The easiest way to compute this Jacobian is to write our C matrix as a GL(A) 
transformation acting on an 'old' matrix 



C 



(30) 



/ 



G 



old 



/ * 

* 

* 
* 



1 









1 









1 








1 / 



(31) 



so that 



C. 



new 
aa 



jf(C new )G 



'Old 

/3a 



(32) 



We know that the measure is simply ^ fc (™- fc )(7 old ) so we can compute the Jacobian in terms of 
the new variables using J _1 C new . It is straightforward to compute this Jacobian in general, 
which we have done in the appendix. In our case, the Jacobian is 



(C4i — C2iC42)(c63 — C^Cq^)(c^ — C§<jC&§)(c27 — C^C2%) 



(33) 
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The product of 4 x 4 determinants in the denominator of the integrand of £ 8 ,4 is 

8 

J J A = (Zl(c 6 3 - C43Q34 + X 3 C 42 Cu))(Cs5(c63 ~ C 43 C 6 4 + X 3 C A2 Cqa)) ■ ■■ 
i=l 

= (C 4 1 - C2lC 42 ) 2 (c 6 3 - C4 3 C 64 ) 2 (C 8 5 - C 65 C 86 ) 2 (c 27 - C 87 C 28 ) 2 

xc4ic 63 c 85 C27 • xix 3 x 5 x 7 + 0(x 5 ) (34) 

Note that the four factors on the first line are squared, but one of each will be canceled by 
the Jacobian. Taking this into account, we see that £ 8)4 takes the simple form 

dxidx 3 dx 5 dx 7 



8,4 



x x x 3 x b x 7 



x f d 12 c u 5 A \\C aa W a ) 

J C4lC63 c 85 c 27( c 41 — C2lC42)(c63 — C 43 C^){c^ — Cq^,Csq)(c27 — C 8 7C2 8 ) 

The contour integral over the x immediately sets them all to zero, so we have neglected 
higher order terms in these variables. The denominator is precisely what we get from the 
denominators of the four MHV amplitudes at the corners of the one-loop leading singularity 
(i.e. the 'box coefficient'; note that four c parameters have been eliminated using r variables). 

Now we can Fourier transform from twistor space back to momentum space. The most 
general way to do this is to write 

J d 2 p a e^^5\C aa W a ) = 5 2 (C aa \ a ) J d 2 Pa 5 2 (\ a - C aaPa ) (36) 

so now the c variables must satisfy 

C aa \ a = and A a - C aa p a = (37) 

where the p a are auxiliary spinor variables. Clearly the first set of equations is linear in the 
C variables. However, because the auxiliary p a are free, the second set of equations is in 
general quadratic. Something interesting has occurred, as the entire kinematic structure of 
the leading singularity is encoded in these simple quadratic equations! Note also that any 
multiplicity of solutions will come entirely from these momentum space equations. We expect 
that in general the multiplicity will exactly match the multiplicity of solutions to the 4L cut 
conditions at L loops. 

The procedure that we have followed generalizes to the computation of any one-loop 
leading singularity, with one crucial caveat - in general, the contour integral over the x 
variables will not be so simple. We will generically have a large C matrix, the number of x 
variables will be much larger than the number of external particles, and the residue at x = 
will be highly composite. However, we have a very definite expectation, namely that this 
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residue must equal the product of the four Cmfc denominators. In the appendix we argue 
for the existence of the residue, but we do not know how to compute it and prove that our 
expectation is correct. 



3.3 Back to BCF 




Figure 2: The BCFW Recursion Relations in maximally supersymmetric theories. 



The BCFW recursion relations p2]-[2l] are an extremely efficient method for computing 
tree level scattering amplitudes in a variety of theories. Some key features of these recursion 
relations are that they compute scattering amplitudes using purely on-shell information, and 
that they assemble local amplitudes from non-local pieces. As an example, the 6-pt amplitude 
in Yang-Mills theory is 



^bcfV" = (l + r 2 + r 4 ) 



(46) 4 [13] / 



x- 



L[12][23](45)(56)(p 4 + P5+P 6 ) 2 

(38) 



(6|5 + 4|3](4|5 + 6|1]_ 

when computed with BCFW (where r cyclicly permutes the external particles i — > i + 1). 
Note that the factor on the second line has unphysical poles, and therefore it could never 
come from the Feynman diagrams of a local theory. One of the main motivations underlying 
the discovery of L n ^ was to find a way to explain how local amplitudes arise from non-local 
pieces. 

The BCFW recursion relations were originally discovered by Britto, Cachazo, and Feng 
[TT] in a study of the IR equations as applied to one- loop leading singularities [52]- [53]. This 
means that each term in the recursion relations is a one-loop leading singularity, so we can 
use our techniques from the previous subsection to identify the contours of integration in 
C Ut k that correspond to tree amplitudes. 

To be more specific, we want to look at one-loop leading singularities with (ni, ki) = (3, 2) 
and (ri2, /C2) = (3, 1), or in other words we take these two neighboring corners of the box to 
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be MHV and anti-MHV 3-pt amplitudes 




where particles 1 and 2 correspond to the analytically continued particles in the BCFW 
recursion relations, and the small unlabeled circles represent projective W variables to be 
integrated over. Using our solution from the previous subsection, we find a C matrix in Cn,k 
of the form 



C 



( Cll! 








V o 



1 

C21 2 















C 2,2 R 
C a,2 R 

C a,2 R 
C I,2 R 








C 2,(n-l) fl 
C a ,(n-1) R 



C a ,2 L 






C a,(n-1) R 
Cl,(n-l) R C I,2l 
C a,2 L 



Ca,(n-1) L \ 






Cl,(n-1) L 
c a,(n~l) L 







(39) 



c a,2 L ■ ■ ■ Ca,(n-1) L / 



where we have indexed most of the c's with L and R to show that these belong to the usual 
M L and Mr of BCFW, and we have used a label / for 'intermediate' for the one overlapping 
row. It is worth noting that this matrix structure is not so surprising - it is perhaps the first 
thing one might guess. The BCFW form of the amplitude is being represented by two blocks 
that correspond to M.^ and M.r and which share a row that corresponds to the intermediate 
particle. 

Let us check our general formula with a few examples. The simplest example is the 
computation of an MHV amplitude by BCFW; for this case we would find a C matrix 



C = 



* 1 * 
* * * 



(40) 



None of the sub-determinants (/, /+ 1) vanish, which is exactly what we would expect for the 
C matrix of an MHV amplitude. The GL(2) symmetry has not been fully fixed because we 



17 



have yet to use the tl projectivity parameter, we can use it to obtain the fully fixed matrix 

\^ o * i * ■■• * y 

from which one could compute the MHV amplitude in momentum space. 

As another example, consider the 6-pt NMHV amplitude. One of the terms used to 
construct it comes from applying BCFW where M.l and M.r are both 4-pt amplitudes. In 
this case we would take the C matrix to be 

/*100**\ 
C= 0***00 (42) 

loo****/ 

We see that only one of the determinants (1,1 + 1,1 + 2) vanishes, namely the one with 
1 = 5. This is precisely what was found in pQ. As another example, consider again the 6-pt 
amplitude constructed from a 5-pt and a 3-pt MHV amplitude, this would have C matrix 

/ * 1 * \ 

C= 0****0 (43) 
\00****y 

so we see that the 1 = 3 determinant vanishes. 

Using these results one can recursively identify the contours of integration that correspond 
to tree amplitudes. A C matrix of the form that we have identified in this section will give 
terms that can contribute to tree level amplitudes as long as the contours of integration for 
cl and cr are chosen to give components of tree level amplitdues. To obtain the full BCFW 
recursion relations one simply sums over the sets L and R with appropriate contours for the 
sub- Gr assmannians . 

The analysis of this subsection and the last is one-half of a constructive proof that all 
one-loop leading singularities are contained in C n> k for some choice of contour. This follows 
because (1) we have (recursively) shown that all tree amplitudes are contour integrals via BCF 
and the IR equations, and (2) we have identified the contour for one-loop leading singularities 
given the sub-contours for the four tree amplitudes at the corners of the box. Our method 
also partially explains why terms that come from non-adjacent BCFW deformations do not 
arise as residues of C n ^ - due to color ordering, these terms cannot be written as one-loop 
leading singularities and so they cannot be found among the residues. However, our results 
are only half of a proof because we have not explicitly computed the residues themselves, but 
only shown how to obtain the appropriate contours of integration. 

3.4 Higher Loops and General Patterns 

In the previous sections we identified the contours of integration for all tree amplitudes and 
all one-loop leading singularities. These results immediately apply to an infinite class of 
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Figure 3: A diagram representing the infinite class of leading singularities that can be iden- 
tified by applying our one-loop analysis recursively, expanding the objects at the corners of a 
one-loop box into new boxes. The heavy black dots are four point amplitudes linking boxes 
together at their corners. 



leading singularities - those that can be constructed by attaching 'boxes' together at their 
corners. This follows because we can recursively interpret each of the 4 sub-matrices of Cn,k 
as one-loop leading singularities themselves. Thus the general statement is that this type 
of leading singularity at L loops corresponds to a configuration where Cn,k is broken up 
into 3L + 1 submatrices following the pattern of equation fl22l) . This structure of leading 
singularity can be visualized as in figure [3] where the black dots represent the joined corners 
and the little tick marks at the other corners represent external particles. 

These sorts of leading singularities always correspond to block diagonal C matrices (note 
that the fact that the blocks lie on the diagonal is itself meaningless because we are free to 
cyclicly translate all of the columns) . This makes sense based on the topological structure of 
the loop diagram, because beginning at any point on the diagram one can follow propagators 
and "walk" from tree amplitude to tree amplitude, encountering every propagator and tree 
amplitude in cyclic order. For more general topologies this would not be possible - one would 
inevitably miss some tree amplitudes and propagators. 

We can write more general leading singularities in twistor space using the method of 
section 2. As a first example we can consider the diagram of figure HI Computing this 
diagram in twistor space is straightforward, since again we only need to integrate over delta 
functions. We will not go through the computation in detail or consider the possible subtleties 
that can arise when the various tree amplitudes at the corners have too few delta functions 
(i.e. for very small and hi). We will only give the generic result because our goal is to 
explicate the pattern of how leading singularities correspond to various sub-structures in the 
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k, n Grassmannian. 




Figure 4: This diagram shows a 2-loop leading singularity and the corresponding points in 
the Grassmannian to which it corresponds. The rectangles in the pictured matrix correspond 
to its non-zero entries, and the adjacent boxes share a single row. 

The easiest way to compute this leading singularity is to first take account of the propa- 
gators around the borders of the box and pentagon and only then integrate over the single 
W variable corresponding to the propagator shared between the box and the pentagon. The 
first step gives a structure in Cn,k that is block diagonal as in the one-loop case except with 
7 blocks instead of 4. The second step eliminates a row and column, with the result that 
two non- adjacent blocks now share a row. This can be pictured as in figure HI where we 
have explicitly displayed the C matrix structure that arises when this leading singularity is 
embedded in Cn,k (the regions outside the boxes are filled with zeroes). This analysis can be 
generalized to another infinite class of leading singularities made up of boxes and pentagons 
that are chained together along various sides in such a way that there are 4L propagators at 
L loops. 

Although we have given the general case above, we have also explicitly checked our results 
for the case n = 12, k — 6, ie for an N 4 MHV amplitude. In that case the full C matrix for 
the 2-loop leading singularity 




5 H 



20 



takes the form 



C 



/ 


Cl,l 


1 


























Cl,ll 


Cl,12 


\ 




a 




c 2,3 


C2,4 


c 2,5 




































c 3,4 


c 3,5 


c 3,6 


1 




































c 4,6 


c 4,7 


c 4,8 


1 















a 


b 


1 














C5,8 


C5,9 


C5,10 










\ 


a 


b 


1 




















C6,10 


C6,H 


C6,12 


/ 



(44) 



This two-loop leading singularity has a very complicated kinematic structure in momentum 
space. By this we mean that when one solves the 8 quadratic equations that force the 8 inter- 
mediate propagators on-shell, the solution involves elaborate double square roots of kinematic 
invariants. When £ 12 6 is transformed to momentum space one obtains the equations 



C aa ^a — and A a — C aa p c 







(45) 



for the Cij and a and b variables, where p a are auxiliary spinors that must be solved for and 
eliminated. We have checked explicitly^ that with our C matrix structure these equations 
give precisely the kinematic structure of the leading singularity. This is an extremely non- 
trivial check of our methods and of the claim that the residues of L n ^ are in fact leading 
singularities. 



f 




10 



\ 



li 



12 



Figure 5: An example of a 4-loop leading singularity and the associated subset of the Grass- 
mannian. The two blue rows are identical, while the two red rows are identical up to an 
overall factor each. 



More interesting cases arise at 3-loops and beyond where we have the possibility of tree 
amplitudes that are entirely internal to the loop diagram. An example of this phenomenon 
is given in figure [5j Here again we have computed the kinematics of the object in twistor 
space by first accounting for the propagators along the boundary and then integrating over 

5 with the help of Jacob Bourjaily 
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the W variables that link the D along the boundary with the internal tree amplitudes. In 
the C matrix structure pictured in figure El the two blue rows are identical, while the red 
rows are identical up to an overall factor each. We see again that the topology of the loop 
graph is reflected in the structure of the subspace of the Grassmannian. 

There are many possibilities for further exploration here, and it may even be possible 
to categorize and understand all of the possible topologies. Other natural goals include 
understanding in a more concrete way how the twistor space structure gives rise to the 
appropriate momentum space kinematics, and understanding whether all residues of C n ^ are 
leading singularities at all loops. It is exciting to note that for any leading singularity our 
methods will give some sub-matrix structure within Cn,k- Thus with one class of exceptions, 
we have implicitly shown that all leading singularities arise from Grassmannian kinematics. 

The exceptions are the so-called "composite leading singularities" [2], which seem to be 
important in obtaining the full loop amplitudes. These are diagrams at L loops with fewer 
than 4L explicit propagators which nevertheless give rise to leading singularities. The classic 
example is the diagram 




where we have shown the series of cuts and manipulations that one can perform in order to 
obtain the leading singularity. The naive translation of this diagram into twistor space would 
seem not to give rise to a leading singularity, but to a product of tree amplitudes integrated 
over one free variable. Clarifying the role that these sorts of leading singularities play in 
constructing general loop amplitudes is an important goal for future work. 

4 Conclusions and Future Diretions 

We have shown how any given leading singularity of the M = 4 SYM S-Matrix can be 
identified among the residues of the Grassmannian contour integral Moreover, we have 
seen that there is a simple and physical pattern for how the various leading singularities 
appear, so that the structure of the perturbation series is reflected in various subspaces 
within the Grassmannian. Let us now consider some directions for future work. 

• Evaluating the Residues The only piece missing from our argument is a method for 
calculating the residues themselves in general - if this could be established, then our 
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argument would become a proof that all leading singularities are residues of £ n ,&- Not 
only is this a precise mathematical problem, but we know the answer ahead of time 
- for instance, at one- loop we know that there must be a residue of £n,k containing 
four smaller C l n . k ., and we know that the actual value of the residue is given by the 
products of the denominator factors from the four smaller L l n . k .. However, this question 
remains both non-trivial and interesting, as the determinant factors that make up the 
denominator of L n ^ make up the 'Grassmannian Dynamics'. Also, the computation of 
multi-variable composite residues is in general a difficult mathematical problem |54j - 
[56], so we expect that the special form of the denominator must play a crucial role. 
It will also be interesting to understand the converse statement, that all residues are 
in fact leading singularities, and perhaps to reverse our logic and formulate a recursive 
'derivation' of C n ,k- 

• Composite Leading Singularities As discussed in section 3.5, we do not have a 
twistor space picture for the composite leading singularites, which arise from diagrams 
at L loops that have fewer than 4L explicit propagators. These diagrams seem to play a 
role in the construction of the full S-Matrix [2], so it may be important to identify them. 
Another possibility is that they are somehow always associated with, or algebraically 
identical to, the more natural leading singularities that we have already identified. This 
is a pressing issue if we hope to unite C n ± with the actual loop integrals to construct 
the full S-Matrix of the M = 4 theor}fl. 

• Kinematic Structures In our two loop 12-pt N 4 MHV example we saw how a par- 
ticular subspace within the Grassmannian automatically encoded the solutions to the 
8 quadratic equations that arise when we 'cut' 8 loop propagators and force them on- 
shell - a rather non-trivial feat. It would be interesting to systematically understand 
how very complicated momentum space kinematics can be encoded by Grassmannian 
subspaces. This may be of particular interest because the 'Grassmannian Kinematics' 
may generalize beyond the M = 4 theory even if the 'Grassmannian Dynamics' (the 
denominator structure and specific residues) does not. 

• Residue Theorems Now that it is possible to identify leading singularities within 
C ni k for very general n and k it will be interesting to try to study the appropriate 
residue theorems [51]- [56] in a systematic way. As we saw in [1], we expect that these 
residue theorems encode the locality of the S-Matrix by enforcing that scattering am- 
plitudes only have physical poles and obey the IR equations. It would be interesting 

6 Unless of course there is some direct, once-and-for-all solution to this problem, as the existence of the 
Wilson Loop/ Amplitude correspondence and dual conformal invariance [22]-[3S] might be taken to suggest. 
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to understand these facts in greater generality and at higher loops. 

• Yangian Symmetry Although the dual conformal invariance |32|-[45| of C n ^ has 

been shown in [30], [31] . an additional miracle occurred, namely that L n ^ was found 
to be proportional to £ n ,fc-2 written in a "momentum twistor space" [50]. In order to 
better understand this miracle, and also because the Yangian generators [41] become 
extremely natural in twistor space, it would be interesting to directly understand the 
Yangian symmetry of This is not so easy because it is only the residues of C n ^ 
that are Yangian invariant; the integrand itself certainly is not. 

Many of the ideas in this paper were inspired by the 'Hodges diagrams' of [TT]-|14]; 
we have made minimal use of them mostly because they would be unfamiliar to most 
readers. Previously, Hodges diagrams have only been used to represent tree amplitudes, 
but our method of writing leading singularities in twistor space shows that one could 
equally well use Hodges diagrams to represent loop-level information (in fact Hodges 
diagrams can enumerate all leading singularities). It has been shown [57] that the 
Yangian symmetry of scattering amplitudes in the Af = 4 theory can be seen via a 
simple induction argument applied to Hodges diagrams. It would be interesting to try 
to extend this argument to all leading singularities. 

• Non-Supersymmetric Theories At one-loop, scattering amplitudes in theories such 
as pure Yang-Mills cannot be characterized by their leading singularities, but require 
the specification of so-called triangle and bubble coefficients and also rational terms 
that have no 4-dimensional unitarity cuts. In 

The Hodges diagram techniques of [H]-[ll] are equally applicable to tree amplitudes in 
Yang- Mills theories without supersymmetry. Using the methods of section 2, it should 
be possible to write triple and double cuts in twistor space, and perhaps with a bit of 
cleverness one could isolate the actual triangle and bubble coefficients. Experience has 
shown that scattering amplitudes come back from twistor space in new and improved 
forms, so it might be useful to attempt to compute pure Yang-Mills amplitudes in this 
way. 

• Building Full Amplitudes It seems reasonable to interpret the very existence of 
£ nt k as an indication of the importance of leading singularities, so it is very important 
to understand if there is some simple way of computing the actual S-Matrix from its 
leading singularities beyond one-loop. 

Another very exciting direction would involve combining the Wilson Loop, which has 
been conjectured to compute MHV amplitudes to all orders [M]-[37J, and which 
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in the form of [30] actually builds non-MHV amplitudes from MHV amplitudes us- 
ing Momentum (or dual conformal) Twistors [50j . More generally, techniques from 
integrability [59], [60] may shed light on 

• Implications for Gravity? A holy grail and initial motivation for much recent work 
has been the hope of finding something like L n ^ for Af = 8 Supergravity [61]- [63], a 
theory whose perturbative S-Matrix may also be determined by its leading singularities 
[3], and may be finite [51]- [Slj. If found, such an object could be viewed as a holographic 
description of flat spacetime. 

The pattern of leading singularities within C n> k gives us hints for how something like 
£ n) fc might work for Af = 8 Supergravity Leading singularities seem to be equally 
important in Af = 8 as in Af = 4 [3], so if a direct analogue of exists for gravity, 
we might expect it to have the same sort of topological and recursive structure as 
we have found for the Af = 4 theory, except without color ordering. It may make 
sense to ask questions along the lines of "does there exist a manifold containing the 
gluing of four gravitational tree amplitudes in all possible permutations?". Also, we 
know from its non-conformal nature, from the fact that the gravitational 'charge' is 
energy-momentum, and from explicit checks that leading singularities in Af = 8 cannot 
be characterized with as few kinematical variables as those of Af = 4, and this again 
points to a some new and different space for a dual description of gravity. 

Note Added: During preparation of a companion paper to this work, an interesting new 
paper |S5]/ appeared which has some overlap with this work. 
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A The Residues of Cn,k 

A.l Jacobians 

Recall that 



d kxn C n 



(12 - - - fc) (23- •• (£;+!) 



(nl---(k-l) 



l[5^(C aa W a ) 



(46) 



a=l 



is invariant under GL(k) transformations that take C aa — > L£Cp a . This is a redundancy 
of description, analogous to the gauge symmetries necessary to provide local descriptions of 
massless spin 1 and spin 2 particles (in our case the redundancy makes the cyclic permutation 
symmetry manifest). This redundancy must be eliminated before we can compute leading 
singularities. 

We have introduced a new gauge fixing for this GL(k) redundancy, so in this section we 
will compute the relevant Jacobian. Perhaps the most difficult issue is coming up with a clear 
notation for these large matrices, so we will refer throughout to an example in the hopes that 
the general case is clear. 

With the 'canonical' gauge fixing of pQ, where the C matrix is fixed so that some k of 
its columns form the k x k identity matrix, the Jacobian is 1. Since our gauge fixing is very 
similar to this one, it will be easiest to compute our Jacobian by transforming from this 
gauge fixing to our own. 

As a rather general example to keep in mind, a C matrix with the 'old' gauge fixing would 

be 



a 



old 



/ * 

* 
* 



X X 
X X 
\ X X 



1 

1 











x 

X 

* 
* 

X 
X 





1 






X X 
X X 
X X 









1 
1 




0\ 




x 






1 1 



(47) 



whereas with our gauge fixing we will take 



( Cn C12 1 X15 Xiy 

C21 C 22 C 23 1 X 2 5 X 2 7 X 2 8 

C31 C 32 C33 C35 1 X37 X 38 

X41 C45 C 46 C47 C 48 








1 



Ci,!! 01,12 \ 





£51 x 52 
\ x 6 i x 62 



c 57 c 58 c 59 





£4,11 

1 £5,11 



c 67 c 68 c 69 06,11 






1 



(4£ 
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It is easy to find the GL(k) transformation that relates these two matrices. We simply take 



J(C new ) c f with J(C new ) 





1 














Cl,ll 






C23 


1 


















C 3 3 





1 





















C46 


1 





















c 59 


1 







V 











C69 





1 


J 



(49) 



We have emphasized that J is a function of the C ncw variables, so that J -1 • C Dew also 
depends entirely on these variables. Now we can compute the Jacobian from the equation 
j-i . c ncw = C old . Taking d of both sides and then multiplying by J gives 



old 



(50) 



Since C n ,k is invariant under global GL{k) transformations, this last multiplication with J 
drops out of the overall Jacobian, which we can now compute directly from the left hand side 
of the equation above. It is amusing that this equation makes it manifest that J is a GL{k) 
"gauge field". Also note that this equation is completely general, and does not depend on 
any of the details of our particular illustrative example. 
Now the measure comes from taking the wedge product 



(51) 



and the variables that do not appear in J can be factored out. This means that the only 
columns (values of a) that produce a non-trivial Jacobian are those where there are extra 0s 
in C ncw . Thus the Jacobian is 



K 



1=1 



dJ~ 



dd 



(52) 



where % labels the K entries in C ncw that have been set to zero by our gauge fixing, and at 
and ctj are the corresponding columns and rows. This formula simply reduces to a product 
of minors to various powers; in the case of our example the Jacobian is 



(C32 - C12C33XC45 - C35C 46 ) 2 (c 68 - C4 8 C 69 )(ci i ii - C^nCi^)' 



(53) 



In general, with our specific gauge fixing, the Jacobian is given by a product of four (ki — 
1) x (hi — 1) minor determinants, each raised to the power — 1. These are the right-most 
minors in each of the (rij — 2) x (hi — 1) sub-blocks corresponding to the four corners of the 
one-loop leading singularity (box), as can be verified by a straightforward computation. 
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A. 2 Existence of Tree and One-Loop Residues 

In this appendix we will argue for the existence of the residues of Cn,k that give rise to 
the block structure of the C matrix corresponding to the four D . We will refer to the 
Grassmannian coordinates that we wish to set to zero as x variables, as pictured in the C 
matrix of equation (|48|) . We will show that the denominator of C^,k vanishes to high enough 
order in the x variables for the point x = to be a residue. 

To begin let us count the number of x variables, noting for convenience that N+8 = ^ Hi 
and K + 4 = £\ fcj. There are NK — ^ k^rii — 2) entries in C outside of the sub-matrices 
corresponding to the C\ but K + Yl(h ~ ~~ are se ^ to zero once we fix the GL(K) 
redundancy, so there are 

4 

N X = (N- K)K + 4 - Mm - k i)l ( 54 ) 

i=l 

x variables in total. Now we need to show that 

D = (12...K)(23...K + 1)...(M2.. Jif - 1) (55) 

has no terms of lower order lower than this in the x variables. Another way of saying this is 
that we want to prove that the denominator, considered as a polynomial in the x variables, 
is to leading order homogeneous and of degree N x . 

It suffices to examine how the rank of the K x K matrices appearing in D depends on 
the x variables. Specifically, we would like to consider how the sum of the ranks of these iV 
matrices changes when x take generic values versus when all x — 0, since this tells us the 
order of D as a polynomial in the x. For each % = 1, 2, 3, 4 there are k{ — 2 rows full of xs (or 
zeroes) that are each of length iV + 2 — rii and also four rows of length N + 4 — rii — n^+i. 
The presence of each row increases the order of D in the x variables by the length of the row 
minus K — 1. However, there is an additional effect near the corners of the D sub-matrices 
because a linear dependence in either the rows or the columns of a matrix will decrease its 
rank. This contributes (hi — l)(fcj — 2)/2 at two corners of each of the four sub-matrices, 
giving a total 

4 

{N - K + 3)K + 8 + i( k i -l){ki-2)- hm] (56) 

i=\ 

This is precisely equal to the number of x variables N x that we counted above. Without a 
better understanding of the precise definition of the residue we cannot conclude that it exists, 
but our argument makes it very plausible. 
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B All NMHV Residues 



Now we will give a solution for all the residues of C n ^. By a solution we mean an explicit 
identification of every residue of the contour integral 

This is a multi-dimensional contour integral over a G(3,n) Grassmannian; it is useful to 
count the number of integration variables in order to see the best way to label the residues. 
After eliminating the GL(3) redundancy of the Grassmannian, £ ra3 becomes an integral over 
3n — 9 variables. When we Fourier transform from twistor space back to momentum space, 
we produce 2n delta functions, but 4 of these turn into the momentum conservation delta 
function, so there are only 2n — 4 independent constraints. After these constraints have been 
taken into account C n $ reduces to a contour integral over (3n — 9) — (2n — 4) = n — 5 free 
variables. The denominator of the integrand is simply a product of n 3 x 3 determinants, and 
on the delta function constraints these are each linear functions of the n — 5 free variables. 
Thus a single residue can be specified by listing the 5 determinants that are not set to zero 
at the residue of the contour of integration. 

It is easiest to think of the solution as being given by this diagram 




Those familiar with jT4] may note that this is a 'Hodges diagram', but knowledge of these 
diagrams is not essential to understand what follows. The particles are labelled by an integer 
from 1 to n, and A, B, C, D, E can be any increasing set of integers in this range. We are 
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representing these particles at the vertices of the pentagon with Z = (A, jl, fj) twistors, while 
all of the other particles, which are not explicitly drawn, are most naturally taken to be 
W = (A, fi, i]) twistors. This is simply a choice of basis and is not physically meaningful, 
but it will be useful in what follows. What the diagram means is that we take the anti- 
MHV 5-pt amplitude M 5 (Z A , Z B , Z c , Z D , Z E ) and up to five MHV amplitudes such as 
A4(Za, Wa+i, Zb) and simply multiply them. The number of particles in each MHV 
amplitude is fixed by differences such as B — A; MB = A + 1 then there is no MHV 
amplitude on the AB side of the pentagon. 

What does the diagram mean physically? It turns out that this is the most general 
object that one can get from applying the BCFW recursion relations to compute NMHV 
tree amplitudes. Our claim is that this diagram is precisely the residue that we would label 
{A, B, C, D, E}, where eg A represents the determinant (A — 1, A, A + 1). Let us now show 
this explicitly. 

First we will write £„ i3 in a basis where particles A, B, C, D, E are represented by Z and 
the others are represented by W in order to facilitate comparison with the diagram. We fix 
the GL(3) redundancy of C n ,z by setting columns A, B, and C to the identity matrix. Next 
we Fourier transform these particles to the Z basis, giving 



-n;3 



J ?T7 ^ " 1 ; ^ 7 —^ e iCijWi - Zj (58) 



d 3n - 9 c u 

(123) •••(«- + l)---(nl2)' 



where J = A, B, C. Now we can Fourier transform particles D and E to the Z basis as well, 
giving 

£ W (123, —i X'Z i) ■ ■ • m e "" w " £ ' sl,l(ZD + °°> z >w z * + ^ < 59 » 

The residue of interest is obtained by setting (I — 1, /, 1+ 1) = for all I ^ A,B, C, D, E. We 
will now see that the diagram can be written as an integral over the same qj variables with 
the same structure of delta functions as £ n ,3- The vanishing of the claimed determinants will 
be guaranteed by the structure of the diagram. 

The central pentagon of the diagram is simply an anti-MHV 5-pt amplitude. In accord 
with our choice of variables for £ n 3 let us represent it in the all Z basis with its C aa matrix 
fixed to be 

c p = ( c ) a c r f c I ? ) (60) 

V C EA C EB C EC u 1 J 

where the p index indicates that these are the c's in the pentagon. Now we can write the 
pentagon as 

(AB)(BC)...(EA) S ^ Z ° + 4u*>)*»lZ* + 4**) («> 
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and the delta functions match up with our representation of £ n; 3- 

Our diagram represents the product of this pentagon with the five MHV amplitudes that 
are attached to its edges. We can write each of these amplitudes as a copy of C m -,2- For 
example, the MHV amplitude attached to A and B can be written with a C AB matrix 

qAB _ f 1 CA+1,A Ca+2,A ■■■ C B -1,A \ 
\ Ca+1,B Ca+2,B ■■■ Cb-1,B 1 / 

so that the MHV amplitude itself takes the form 

/AriAB 
(AA+l)...(B-l,B)(B,Af [b6) 

where J = A,B and i runs from A+l to B—l. Similar expressions obtain for A4 BC , A4 CD , A4 DE , 
and M EA . 

Now we can see that with our choice of "gauge fixing" of the various GL(2) and GL(3) 
Grassmannian redundancies of description, (/ — 1, /, / + 1) = for A < I < B but that 
(B — 1,B,B + 1) does not vanish. The former statement follows from the fact that ci t c 
does not exist, so in other words ci t c = by definition. Since the determinant factors are 
linear in cj c, they vanish. The latter statement follows by direct evaluation - (B — 1, B, B + 
1) = cb~i,aCb+i,c which can be seen to be non-vanishing in momentum space by a direct 
computation. 

Both the diagram and £ n>3 are independent of the choice of Z or W basis and the "gauge 
fixing" of the various GL{2) and GL(3) redundancies. With different gauge fixings it would 
be clear that the determinant (/ — 1, 1, 1 + 1) = for all I ^ A,B, C, D, E. Since we have 
made no assumptions that break the symmetry between A, B, C, D, E except for the choice 
of basis and "gauge" , we can conclude that the diagram corresponds to the claimed residue 

of £„,3- 
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